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FOREWORD 


McDonnell Douglas Helicopter Company (MDHC) has been conducting a study of finite element 
modeling of helicopter airframes to predict vibration. This work is being performed under U.S. 
Government Contract N'ASl-17498. The contract is monitored by the NASA Langley Research Center, 
Structures Directorate. 

This report summarizes the procedure used at MDHC for predicting coupled rotor/fuselage vibrations 
with an application to the AH-lG two-bladed rotorcraft including comparisons with flight test 
vibrations. Key NASA and MDHC personnel are listed below. 

NASA Langley McDonnell Douglas Helicopter Company 

Janice H. Clark, Contracting Officer M. Toossi, Project Manager 

Joseph VV. Owens. Contract Specialist R. King, Program Manager 

John H. Cline, Technical Representative K. Sangha, Rotor Dynamics 

Raymond G. Kvaternik, Leader *J. Shamie, Rotor Dynamics 

Rotorcraft Structural Dynamics Group 


‘Now with Hughes Aircraft Co. 
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INTRODUCTION 


INTRODUCTION 


The need for developing and implementing improved capabilities for accurately determining 
rotor /airframe response loads is well recognized. Presently, reduction of rotor /airframe loads to an 
acceptable level in many cases required intensive, costly, time and payload consumming design 
modifications. The NASA Langley Research Center is sponsoring a rotorcraft structural dynamics 
program with the overall objective to establish in the United States a superior capability to utilize 
finite element analysis models for calculations to support industrial design of helicopter airframe 
structures. Viewed as a whole, the program is planned to include efforts by NASA, Universities, and 
the U.S. Helicopter Industry. In the initial phase of the program, teams from the major U.S. 
manufacturers of helicopter airframes will apply extant finite element analysis methods to calculate 
static internal loads and vibrations of helicopter airframes of both metal and composite construction, 
conduct laboratory measurements of the structural behavior of these airframe, and perform 
correlations between analysis and measurements to build up a basis upon which to evaluate the results 
of the applications. To maintain the necessary scientific observation and control, emphasis throughout 
these activities will be on advance planning, documentation of methods and procedures, an thorough 
discussion of results and experiences, all with industry-wide critique to allow maximum technology 
transfer between companies. The finite element models formed in this phase will then serve as the 
basis for the development, application and evalutation of both improved modeling techniques and 
advanced analytical and computational techniques, all aimed at strengthening and enhancing the 
technology base which supports industrial design of helicopter airframe structures. Here again, 
procedures for mutual critique have been established, and these procedures call for a thorough 
discussion among the program participants of each method prior to the applications and of the results 
and experiences after the applications. 

This report is a description of an analysis method, plan and results obtained in calculating coupled 
rotor, airframe flight vibration levels for a helicopter in steady state forward flight The method used 
is MDHC’s Rotor Airframe Comprehensive Aeroelastic Program (RACAP). 
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PHILOSOPHY OF ANALYTICAL APPROACH 


RACAP is a comprehensive rotorcraft design and analysis tool which has been formulated with 
state-of-the-art philosophy in helicopter aeroelasticity. In the structural representation of blade 
elements emphasis is placed on completeness of the model. Coupling of the rotor with the fuselage is 
performed in an accurate manner. The aerodynamic model has been developed with provision for 
updating airloads, unsteady aerodynamics, and inflow modelling. A modular approach is used for ease 
of interchangeability and improvement of analytical methods. 



PHILOSOPHY OF ANALYTICAL APPROACH 



• “COMPLETENESS” IN STRUCTURAL 
REPRESENTATION OF BLADE ELEMENTS 

• PROVISION FOR ACCURATE COUPLING OF 
ROTOR WITH FUSELAGE 

. AERODYNAMIC MODULE WITH PROVISION FOR 
UPDATE OF 

• AIRLOADS 

• UNSTEADY AERODYNAMICS 

• INFLOW MODEL 

• MODULAR APPROACH FOR EASE OF 
INTERCHANGEABILITY AND UPGRADE OF 
ANALYTICAL METHODS 
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3.0 OVERALL CAPABILITIES OF RACAP 



OVERALL CAPABILITIES OF RACAP 


RACAP serves as a comprehensive aeroelastic computational program for prediction of rotor /airframe 
characteristics. Included in the existing capabilities of RACAP are: 1) rotor aeroelastic loads 
prediction, 2) airload calcualtion for use in far-field noise prediction, 3) free vibration analysis, 4) 
fuselage vibration prediction, 6) calculation of performance characteristics, 7) prediction of 
performance and vibration effects of Higher Harmonic Control for all rotor and fuselage configurations 
and 8) bearingless rotor analyses including details of the multiple load paths. 



OVERALL CAPABILITIES OF RACAP 


A COUPLED ROTOR/AIRFRAME FINITE ELEMENT PROGRAM 
FOR: 

• ROTOR AEROELASTIC LOADS PREDICTION 

• AIRLOADS CALCULATION FOR USE IN FAR-FIELD NOISE 
PREDICTION 

• FREE VIBRATION ANALYSIS 

• FUSELAGE VIBRATION PREDICTION 

• PERFORMANCE CHAR ACTERISTICS 

• PERFORMANCE AND VIBRATION EFFECTS OF HIGHER 
HARMONIC CONTROL FOR ALL ROTOR (ARTICULATED, 
TEETERING, HINGELESS, AND BEARINGLESS) AND 
FUSELAGE CONFIGURATIONS 
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4.0 APPROACH 
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HIGHLIGHTS OF RACAP FORMULATION 


HIGHLIGHTS OF RACAP FORMULATION 


Tho approach - rnployed in RACAP performs rotor trim and airloads calculations in the time domain. 
The coupled rotor/airframe analysis is based on coupling a state-of-the-art aeroelastic rotor loads 
modei formulated using the transfer matrix approach, with a NAS I RAN fime element/experimental 
airframe model. Impedance matching at the huh performed with a 6 degroe-of- freedom hub impedance 
matrix in combination with a harmonic balance solution is used in the coupling procedure. The 
structural response of the blades is obtained as a superposition ol the harmonic responses due to 
harmonic loads. Fuselage vibrations are determined by subsequently applying the calculated hub loads 
to the iVASTRAN finite element model. 



HIGHLIGHTS OF RACAP FORMULATION 


RACAP IS BASED ON: 

• ROTOR TRIM AND AIRLOAD CALCULATIONS PERFORMED 
IN THE TIME DOMAIN 

• ROTOR TRANSFER MATRIX APPROACH 

• FUSELAGE 6 D.O.F. HUB IMPEDANCE MATRIX AND 
AIRFRAME DYNAMICS OBTAINED USING A NASTRAN 
FINITE ELEMENT MODEL 

• ROTOR /AIRFRAME COUPLING THROUGH IMPEDANCE 
MATH C TNG AT ROTOR HUB IN COMBINATION WITH A 
HARMONIC BALANCE SOLUTION 

• BLADE STRUCTURAL RESPONSE OBTAINED AS 
SUPERPOSITION OF RESPONSES DUE TO HARMONIC 
LOADS 

• HUB LOADS [SUBSEQUENTLY] USED TO CALCULATE 
FUSELAGE VIBRATION USING FINITE ELEMENT MODEL 
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4.2 RACAP LOAD PREDICTION FLOW 

CHART 
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raoap loads prediction flow chart 

h<> general helicopter aeroelastic problem is <li v id«*<l into six major categories: 

1. Air mass dynamics (induced inflow modeling) 

2. Calculation of aerodynamic loads 

3. Rotor blade dynamics 

1. Blade/fuselage coupling 

■ r >. Determination of control inputs 

b. Fuselage dynamics 


Coupled rotor, fuselage blade dynamic response is obtained in the frequency domain. The harmonics 
are combined into a time history of response that is used in the aerodynamics module to compute new 
estimates of control angles (trim) and aerodynamic loads. The new loads and control angles are used 
m the structural module. The process is repeated iteratively until a converged aeroelastic solution is 

o o 1-ell ii i . 



ROTOR LOADS PREDICTION FLOWCHART 


RACAP 


STRUCTURAL MODULE 
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AEROELASTIC SOLUTION 
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Rotor/Airframe Comprehensive Acroelastic 1* 


rogram 
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TRIM ANALYSIS 


TRIM ANALYSIS 


I he purpose of the * rim analyses is to provide control inputs (mam rotor collective and cyclic pitches) 
and an initial estimate of the rotor airloads for use in RACAP. This is done initially using a rigid blade 
model, in an external program. An isolated rotor trim algorithm is typically used to trim the rotor for 
a given hit and propulsive force using two rotor control variables: the main rotor collective and 
longitudinal cyclic pitches. The trim program uses a lifting line model to evaluate the forces and 
moments acting on the rotor blade. Subsequent to the first aeroelastic iteration, the control settings 
are obtained as a function of the elastic deformation within RACAP, using a Newton-Raphson or 
Secant method approach. A longitudinal thrust and propulsive force trim only is performed. 



'HUM ANALYSIS 


• PROVIDES CONTROL INPUTS AND INITIAL MAIN ROTOR AIRLOADS 

• ISOLATED ROTOR: TRIM ON LIFT AND PROPULSIVE FORCES. 

USES 2 TRIM VARIABLES (MAIN ROTOR 
COLLECTIVE AND LONGITUDINAL CYCLIC 
PITCHES) TO TRIM 

• FORCES AND MOMENTS DUE TO MAIN ROTOR EVALUATED USING 
LIFTING LINE MODEL FOR BLADE AERODYNAMICS 



TRIM ANALYSIS 


anuT^ -h hlade is identical to the other 

Tip-Path-Plane tilt is the primary contributor to trimZd^nlTanl ? he 

solved in the time domain and the Deriodir blade rW , " g B1 ade dynamic equations are 

scheme. The model uses a lifttl Le ^ I i 7 f ° bta ‘ ned Using a "“™ncal integration 

— 

velocity at the rotor disk. The effects of tip-lj l^Zh " M "‘ C * d 

helicopter ZurZ ™" P "T"" ' " rr "' t f ” th ' *«« «d nroments „„ the 

to approach the Wi * kl - *" « tnaooer repaired 

the help of a ounrerLlly evah.aS S m se^iZ“^7T ^ ” Z™"' WM 
periodic blade motion solution is calculated m** * . *, . 1 ie tnmmed st ^te along with a 

provided as initial values for use in RACAP ^ ^ ^ ° a< S ^ harrTlon,cal, y analyzed and are 





TRIM ANALYSIS (continued) 


• MAIN ROTOR MODEL: 

• ISOLATED ROTOR TRIM ANALYSIS USES BLADE FLAPPING 
DEGREE OF FREEDOM ONLY. 

• BLADE DYNAMIC EQUATIONS ARE SOLVED IN THE TIME 
DOMAIN USING NUMERICAL INTEGRATION SCHEME 

• LIFTING LINE (BLADE ELEMENT) THEORY ALONG WITH 
EXPERIMENTAL TWO-DIMENSIONAL AIRFOIL DATA ARE 
USED TO COMPUTE THE AIRLOADS IN THE TIME DOMAIN 

• OPTION TO USE EITHER SIMPLE NONUNIFORM (Glauert’s) 
INFLOW MODEL OR RIGID OR FREE WAKE INDUCED 
INFLOW MODEL 

• ONCE TRIM IS OBTAINED, AIRLOADS ARE HARMONICALLY 
ANALYZED AND ARE PROVIDED AS INITIAL VALUES TO 
RACAP 
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5.2 AIR MASS DYNAMICS (INDUCED 
INFLOW MODELLING) 


INDUCED INFLOW MODELLING 


RACAP has provisions to use three different induced inflow models in its main rotor model. They are 
(i) uniform induced inflow, (ii) simple non-uniform induced inflow (Clauert’s model) and (iii) rigid or 
free vortex wake induced inflow. 



INDUCED INFLOW MODELING 


• RACAP HAS PROVISION TO USE THREE DIFFERENT 
INDUCED INFLOW MODELS 

• UNIFORM INDUCED INFLOW 

• SIMPLE NON-UNIFORM INDUCED INFLOW (GLAUERT’S 
MODEL) 

• RIGID OR FREE WAKE INDUCED INFLOW 
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5.2.1 UNIFORM INDUCED INFLOW MODEL 
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UNIFORM INDUCED INFLOW MODEL 

mom^um^theory!^ ^ ^ mam r ° t<>r d ' Sk ' S 358,11116(1 to be constant and is evaluated using simple 
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UNIFORM INDUCED INFLOW MODEL 


• UNIFORM INDUCED INFLOW MODEL 

• INDUCED VELOCITY OVER THE MAIN ROTOR DISK IS 
ASSUMED TO BE CONSTANT 

• INDUCED VELOCITY EVALUATED USING SIMPLE 
MOMENTUM THEORY 
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5.2.2 SIMPLE NONUNIFORM INDUCED 

INFLOW MODEL 
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SIMPLE NON-UNIFORM INDUCED INFLOW MODEL 

and aft direction of the rotor disk 

around the azimuth. This model incorporates the effect nf y anatlon in the educed velocity 

rolling moments) on the rotor disk. 06 aer °dy na mic moments (pitching and 
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SIMPLE NON-UNIFORM INDUCED INFLOW MODEL 


• INDUCED INFLOW IS ASSUMED TO VARY LINEARLY IN 
THE FORE AND AFT DIRECTION AROUND THE AZIMUTH 


41 




CNJ 



WAKE INDUCED INFLOW MODEL 



VORTEX WAKE INDUCED INFLOW MODEL 

This inflow model uses a deformed ? 

ffwa^ R fl SearCh Laboratories 1). ^ y aXn 6 S Tofflki;f y d ^? loped * the *ASA Division of 
of wake influence coefficients to iteratively generate non f r g C °f dlt,ons > th * model generates a set 

wl°e%tLr^ Wake , Pr0gram has opLTof ^;^:^; ndUC !? Ve, ° CitieS ° Ver the -in 

wake geometry depending on the flight condition Each kl ^ f (deformed) or rigid (prescribed) 

lv°t r Ce rr edShed and trai,ed VOrtic «ty following each Hftt 1S , m ° dele l d 33 a line, with a mesh 

only trailed voticty. The generation of wake geometry t do k'" 6 “ 6 near wake - The wake has 

rotor in a free stream. In a general free-wake modeT tk T ^ * PT ° CeSS simi,ar to start-up of a 
the generation process. Vortex element end oointw’ ^ Wake dements are allowed to freely distort in 

ve,oc,ty , d „ ™ induced velocity - -rir 
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VORTEX WAKE INDUCED INFLOW MODEL 



• USES VORTEX WAKE MODEL ORIGINALLY DEVELOPED BY 


• EACH BLADE IS MODELED AS A LIFTING LINE 

• MESH OF SHED AND TRAILED VORTEX ELEMENTS WITH 

RIGID FINITE CORES (FULL MESH) BEHIND EACH BLADE 
IN THE NEAR WAKE 

• TRAILING VORTEX ELEMENTS (MODIFIED MESH) ONLY IN 

THE FAR WAKE ' 

• WAKE GEOMETRY GENERATED BY A PROCESS SIMILAR 
TO START UP OF A ROTOR IN FREE STREAM 

• OPTION OF USING EITHER A FREE (DEFORMED) OR RIGID 
(PRESCRIBED) WAKE GEOMETRY DEPENDING ON THE 
FLIGHT CONDITION 

• GENERATION OF WAKE INFLUENCE COEFFICIENTS FOR 
ITERATIVELY DETERMINING NON-UNIFORM WAKE 
INDUCED VELOCITIES OVER THE MAIN ROTOR DISK 
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5.2.4 SCHEMATIC OF WAKE GEOMETRY 

ELEMENTS 
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WAKE GEOMETRY MODEL 


each blade (modeled as a HflLtL^^und^ort^x^ followed *by tra ' led VOrtices '^mediately behind 
for the remainder of the wake The vortex ' . d by se * ° f one or more trailing vortices 

strength and core radius along 25 ZtZ “ ^ ^ 
number of revolutions of the rotor is retained for 8 • • De J? en ? 1 “ g on the advance ratio, a chosen 

The induced velocities at the wl wTnta L 17T7 ,h . 6 ' ndU,:ed Velociti “ at the 

Biot-Savart law. P the b ' ade load P° ,nts are evaluated using the 
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WAKE GEOMETRY MODEL 
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GENERAL FEATURES OF THE 
AERODYNAMICS MODEL 


general features of the aerodynamics model 

th r a, T with — <? 

stall, an unsteady aerodynamic model based on thi * f ^ ^ a,r oac * Sa ^ or conditions below 

evaluated using semi-etnpirical Zie “lie f ‘ ^ ‘ he ° ry Dy ™ m>C Sta " <*«*■ « 

flow on the blades, US ' d ‘° “»**• U- o<f«ts of radial 

computations are made in the time domain an d ; , i j t! £ S a " d reVerse flow ’ The air,oad 
The aerodynamics model has the option of using any one^f th^ff f ngid ^ deformations - 

described earlier. If the vortex wake induced infl * Ti - the d,f f erent educed inflow models 
compatible set of induced velocities and blade airloadT Th S US ® d ’ mflow iterations are made to get a 
wth the blade response solution to get a c^^ 
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GENERAL FEATURES OF THE AERODYNAMICS MODEL 


• LIFTING LINE (OR BLADE ELEMENT) THEORY IN 
COMBINATION WITH EXPERIMENTAL TWO-DIMENSIONAL 
AIRFOIL CHARACTERISTICS 

• UNSTEADY AERODYNAMICS BASED ON THIN AIRFOIL 
THEORY 

• SEMI-EMPIRICAL DYNAMIC STALL MODEL 

• SIMPLE FORMULATIONS TO CONSIDER THE EFFECTS OF 
RADIAL FLOW ON THE BLADES, THREE DIMENSIONAL 
EFFECTS NEAR THE TIPS, AND REVERSE FLOW 

• EFFECTS OF RIGID AND ELASTIC BLADE MOTIONS 
CONSIDERED 

• AIRLOAD COMPUTATIONS MADE IN THE TIME DOMAIN 

• FOR THE CASE OF VORTEX WAKE INDUCED INFLOW, 
INFLOW ITERATIONS PERFORMED TO OBTAIN A 
COMPATIBLE SET OF INDUCED VELOCITIES AND BLADE 
AIRLOADS 

• AIRLOAD CALCULATIONS ITERATED WITH BLADE 
RESPONSE SOLUTION 
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UNSTEADY AERODYNAMICS MODEL 



UNSTEADY AERODYNAMICS 

^mpu^ations^shlg^nsteady^hi^ airfoil U th Stea< ^^ a "°f ynimics are tocluaed in the airloads 
along the blades ** r ° t<>rS ' The <*«*= “ f radial «°» 

slope and aerodynamic center are included in th ^ COITectlons for rea ^ A° w effects on the lift curve 
moment. No unsteady Cto am Mud T? ™ U *‘ >y lift lnd 

conditions, provision is ml 5, ** “ 

a^ktSgT^ri “ d^ PitCh T 6 “ Ue d6layed hi « ha ' “*'« o' 

leading edge are included ^ (deptdTng ^ ° f shad fro “ ** 

at Massachusetts Institute of Technology. ? US '” g “ emp,r,cal model originally developed 
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UNSTEADY AERODYNAMICS 


• USE OF THIN AIRFOIL THEORY SUITABLY MODIFIED FOR 
ROTORS IN CONDITIONS BELOW STALL. EFFECTS OF 
RADIAL FLOW ALONG THE BLADES, TIME VARYING 
ONSET FLOW AND CORRECTIONS FOR REAL FLOW 
EFFECTS ON THE LIFT CURVE SLOPE AND AERODYNAMIC 
CENTER ARE INCLUDED IN THE EXPRESSIONS FOR LIFT 
AND PITCHING MOMENT 

• NO UNSTEADY DRAG EFFECTS BELOW STALL 
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6.0 BLADE STRUCTURAL MODEL 
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6.1 DEFINITION OF GLOBAL COORDINATE 

SYSTEM DOFs 
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DEFINITION OF GLOBAL (HUB FIXED) COORDINATE SYSTEM 

rotating) ° r ^ 

blade are defined. These are rotations and dLpl^emLt, in f T P ° S ! tl0n at each station of the 

kinem t d | reCti ° n is elast ic/kinematic torsion w hi le the^x 1 A ^ dkections ' The rotation 

kmema ic. The global coordinate system is defined al X l ? la k dlSp ac * nt » entirely 
v ert cally upward. The constitutive equilibrium al ° ng the blade and Z 

n g ht-han ded system fixed at the shear center of the ^“ patlbdlty ' equations are derived in a local 
global, hub-fixed system in the undeformed geometry. b and then transf °rmed to the 
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DEFINITION OF GLOBAL (HUB FIXED) 
COORDINATE SYSTEM 



FLAP WISE STATE VARIABLES 
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Definition of Global Hub Coordinate System Lag DOFs 



DEFINITION OF GLOBAL (HUB FIXED) 
COORDINATE SYSTEM (continued) 


CHORDWISE STATE VARIABLES 
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definition of 


global coordinate system torsion 


DOFs. 
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DEFINITION OF GLOBAL COORDINATE SYSTEM 
TORSION DOFs (continued) 



67 


00 



PRECEDING PAGE CLANK NOT FILMED 


6.2 GLOBAL-LOCAL COORDINATE 
TRANSFORMATION 
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GLOBAL-LOCAL COORDINATE TRANSFORMATION 
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GLOBAL-LOCAL COORDINATE TRANSFORMATION 


GLOBAL-LOCAL COORDINATE TRANSFORMATION AT STATION i 


L 

• LOCAL COORDINATE SYSTEM IS ATTACHED TO THE 
SHEAR CENTER OF THE DEFORMED BLADE CROSS 
SECTION AT STATION i 

• e T J L ,e F SEQUENCE OF ROTATIONS WITH RESPECT TO 
GLOBAL (HUB-FIXED) COORDINATE SYSTEM 
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6.3 


BLADE ELEMENT REPRESENTATION 
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blade element representation 


displacement and rotation compatibility of a beam element i^Td k"' an<1 moInent e 1»i«britim and 
equations are formulated usintt state variables 7* T r '° Cat f betwee " sta ‘i<™ i and i+1. All 

in the global coordinate system. The (iapwise indThTr ] ' T’ dlSpiacements ’ and rotations) defined 

respect to the neutral axis whereasThe io^d 11 ^‘ Se , ? "T™ 5 are Ca ' Culated 
center. ments are calculated with respect to the shear 
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FLAP WISE EQUILIBRIUM AND 
COMPATABILITY 



flapwise equilibrium and compatibility 



FLAPWI SE MOT TO NT- 

FORCE EQUILIBRIUM: F z = F z + p.z 

MOMENT EQUILI BRUM: Mf = + _ r - 

z DISPLACEMENT- +M /i + Mf + + M F AOx 

•" ROTATION: = f?' _ S it,c9 ' $ " - 2 it , 

«/asfic btitlt—4w 

WHERE FOR STATION ^ 


E/ = 


F z 

r DF 

F z 

r DT 


Mf = 


V1 1F 


Kk 

M AO 


MOUON NAM,C DAMP,NG FORCE DUE TO PITCHING 
INERTIA MOMENT ABOUT THE X-AXIS 

THE X-AXIS 

AERODYNAMIC MOMENT ABOUT THE X-AXIS 

S?HF 1 OFFSEfF 0 RO Y M NAMIC 
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FLAP WISE EQUILIBRIUM AND COMPATIBILITY 
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6.3.2 CHORDWISE EQUILIBRIUM AND 

COMPATIBILITY 
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CHORDWISE EQUILIBRIUM AND COMPATIBILITY 


CHO RD WISE MO TION- 

FORCE EQUILIBRIUM: 
MOMENT EQUILIBRUM: 

X DISPLACEMENT: 

0 L ROTATION: 


^ = Ll F *\ +F ** + F ™ 

Lr ■*', + - «) - F.UXu, - X ( ) 

+Mk+Mh, l+ Mi p + Mi 0l 

S L ‘ t ' S, \L Xi + M ‘"“ - *<♦■*—• M. 
* ~ w «+l ° .+i ~ 9 , +1 


WHERE FOR STATION 


elastic built-in 


F? = 


F x 

r DF 

F x 

t DT 


Mf = 


M/r 


^ao 


A N FR^ntN?^ E A F * r = AERODYNAMIC FORCE 

rafz d,rection AMPING F ° RCE DUE T ° motion [n 

MOTION NAMIC DAMPING F °ROE DUE TO PITCHING 

INERTIA MOMENT ABOUT THE Z-AXIS 

MOMENT OF THE INERTIA FORCES ABOUT THE Z AYK 

DUE TO CROSS SECTIONAL OFFSETS 

AERODYNAMIC MOMENT ABOUT THE Z-AXIS 

MOMENT ABOUT THE Z-AXIS OF THE AERODYNAMIC 

FORCES DUE TO SHEAR CENTER OFFSET 

FROM NEUTRAL-AXIS 
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CHORDWISE EQUILIBRIUM AND COMPATIBILITY 
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6.3.3 TORSION EQUILIBRIUM AND 

COMPATIBILITY 


35 


TORSIONAL EQUILIBRIUM AND COMPATIBILITY 


TO RSIONAL MOTIO N- 


RADIAL FORCE EQUILIBRIUM: 
MOMENT EQUILIBRUM: 

Y DISPLACEMENT: 

9 T ROTATION: 


F? = Fr + F& l+ FL 


ARi 


ARi + Mloi + *Wn 


DTi 


F x 

r i+l 


MT = MT»+MT I + Ml 

(S+i - Z,) - -X,) 

~ ^+' — - Yi+ikintmalie 

»T = >1, - * r ,„ -f», -el 


• +1 

c/asftc built - 


input 


where for i — NHGl-1, 0f nput — input pitch angle (collective, cyclic) 
i * NHGl-1, 0f nput = 0 

and NHGl is the station number giving the location of the torsion hinge 
WHERE FOR STATION t 


F y = INERTIA FORCE , F\ r = AERODYNAMIC FORCE 
MJ = INERTIA MOMENT ABOUT THE Y-AXIS 

MJ f = MOMENT OF THE INERTIA FORCES ABOUT THE Y-AXIS DUE TO CROSS 
SECTION OFFSETS BETWEEN SHEAR CENTER AND NEUTRAL AXIS 
Mar = AERODYNAMIC MOMENT ABOUT THE Y-AXIS 

m ao = MOMENT ABOUT THE Y-AXIS OF THE AERODYNAMIC FORCES DUE TO 
THE CROSS SECTIONAL OFFSETS BETWEEN SHEAR CENTER AND 
AERODYNAMIC CENTER ‘ 

Mqt — MOMENT OF THE AERODYNAMIC DAMPING FORCE DUE TO PITCHING 
MOMENT 


8f> 
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BLADE CROSS SECTION DEFINITION 



BLADE CROSS SECTION DEFINITION 


of variabie p ~ to -*«• ^ 

L ,pa! S 1"^ P °'T f0 ' th r be r WiSe a " d Ch ° rdwise l0Cati0 " ° f “» 1-' center (SC), 

blade root the PA k ^ ° ^ •’ ”® utra * axls ( NA ). and aerodynamic center(AC). At the 

Diaae root the FA is assumed to be coincident with the SC. ; 


2. Arbitrary orientation of local structural principle axes (located at the SC). 

3. Option to input center of rotation (CR) - pitch axis (PA) offset at the blade root (torque offset) 
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6.5 SPANWISE DISTRIBUTION OF 
STRUCTURAL, INERTIAL, AND GEOMETRIC 

PROPERTIES 
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SPAN WISE DISTRIBUTION OF STRUCTURAL, INERTIAL, AND GEOMETRIC 

PROPERTIES 


The transfer matrix method used for the rotor blade dynamic solution requires that the physical 
characteristics of the blade be specified at selected spanwise stations. RACAP implements a 
technique which uses blade distributed data per running inch of span length to determine the 
equivalent blade properties for the model. Equivalent inertia properties are calculated from 

mid-station to mid-station while equivalent structural and geometric properties are calculated 
from station to station. 
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SPANWISE DISTRIBUTION OF STRUCTURAL, INERTIAL, 
AND GEOMETRIC PROPERTIES 


• EFFICIENT AND EASY TO IMPLEMENT PROCEDURE 
FOR GENERATING EQUIVALENT BLADE MODEL 
PROPERTIES FROM GIVEN DISTRIBUTED BLADE DATA 


• EQUIVALENT INERTIA PROPERTIES CALCULATED 
FROM MID-STATION TO MID-STATION 

• EQUIVALENT STRUCTURAL AND GEOMETRIC 
PROPERTIES CALCULATED FROM STATION TO 
STATION 
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7.0 BLADE ROOT END MODELLING 

CAPABILITIES 
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ARTICULATED ROTORS 


P ne, twist, droop, and sweep can be included in the model formulation. 
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ARTICULATED ROTORS 


• COINCIDENT AND NON-COINCIDENT HINGES 

• ARBITRARY SEQUENCING OF ROOT HINGES 

• PROVISION FOR KINEMATIC COUPLING AT HINGE 
LOCATIONS 

• INCLUSION OF GEOMETRIC BLADE PROPERTIES SUCH 
AS PRECONE, DROOP, SWEEP, ETC. 
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TEETERING ROTORS 


Teetering rotors with geometric blade properties such as precone, undersling, built-in-twist, etc. 
can be modeled m a straightforward manner in RACAP. Provisions for including root end 
kinematic couplings such as flap-pitch or pitch-lag coupling are included. A special feature of 
RACAP is that it lends itself well to the treatment of the teetering rotor root boundary 
conditions. Since the solution is formulated in the frequency domain, the root boundary 

conditions are cantilevered for the steady state and even harmonics and pinned for the odd 
harmonics of the blade motion. 
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TEETERING ROTORS 


• PROVISION FOR INCLUDING UNDERSLING, PRECONE, 
DROOP. 

• CONSIDERATION OF ROOT END KINEMATIC COUPLINGS 

• SIMPLIFIED TREATMENT OF BOUNDARY CONDITIONS 
DUE TO CHOICE OF FREQUENCY DOMAIN SOLUTION: 
STEADY STATE AND EVEN HARMONICS MODELED 
WITH CANTILEVERED BOUNDARY, ODD HARMONICS 
MODELED WITH PINNED BOUNDARY 
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HINGELESS/BEARINGLESS ROTORS 


HINGELESS/BEARINGLESS ROTORS 


RACA^ForXZ^r^’T aSSOcia,ed J with «d bearingless rotors can be modeled in 

“ „ C k „ F °Vj 'T leas rotor ’ P reco "«’ droop, and kinematic coupling at the feathering hinge 

whh , h e r r„ r m od n ;, n ; h “ b ***, S “ ffneSS “ d dy " amic charactenstlcs coupied 

analysis^ For bearinglesTrotora th^RA^^^modeUn^rp^mt^tl^multfpL^loat^pat^^^w 

“ b^:: a „Veirb::: P ^ compaiibility at * he a - d — d 
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HINGELESS/BEARINGLESS ROTORS 


• HINGELESS ROTORS 

• GEOMETRIC PROPERTIES SUCH AS PRECONE, DROOP. 

• KINEMATIC COUPLING AT THE FEATHERING HINGE 

• BEARINGLESS ROTORS 

• REDUNDANT LOAD PATHS 

• MODERATE DEFLECTIONS DUE TO BENDING 

• COMPATIBILITY AT INTERFACE 
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BLADE FORCED RESPONSE TRANSFER 
MATRIX SOLUTION TECHNIQUE 
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8.1 BLADE ELEMENT TRANSFER MATRIX 

MODEL 
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BLADE ELEMENT TRANSFER MATRIX MODEL 

discretize the equaUois^f motion^for^he solut^^th 1 ^ flap ‘ ,ag ' torsion structural model is used to 
is classified as ^ ^ dynamics problems - This method 

force quantities are retained as ^ disp . ,acement and 

the aeroelastic response The transfer mat™ t f * tstribution is directly obtained as part of 

state variables defined with respect to the rotating hlled 2^°“ ’ C ° nSiStS * 

Sta H ,i0 " *■ + 1 « •«— by formulating equations 

stations. All aquat^a d"ri“d"n tt hXZd “ Tl' 0 " betw “" ‘ b * <»o 

nonlinear. These equations are linearis a K * co ° rdiI | ate s y ste ni. The resulting equation are 

independent in the rotating system) equilibrium nnsV ^ W " S ° lutlon for tlie stead y (azimuth 
position solution for the direct harmonic nerti h^’ ' 100 and ab ° ut the converged steady equilibrium 
found by a relaxation ^hn7ql whe r e n inc eme^ ^ ^ e ^ uilibriu - P-^ion is 

iteration of the structural response The are derived from the *“t 

iteration solution, and the process repeated until a * mcr j mentaI e< l uat ions is added to the previous 
harmonic solution is obUmadTI S . , f TT* '“'“‘l 0 " “ finally obtained ' Tbe 
about the converged steady solution. ° u ,on ° e armonic incremental equations, linearized 
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BLADE ELEMENT TRANSFER MATRIX MODEL 


• FULLY COUPLED FLAP-LAG-TORSION MODEL 

• HYBRID FORCE-DISPLACEMENT FORMULATION 

• EQUATIONS OF EQUILIBRIUM AND COMPATIBILITY ARE 
SOLVED IN HUB-FIXED GLOBAL FRAME 

• HARMONIC SOLUTION IS ABOUT THE STEADY 
EQUILIBRIUM POSITION 

• DIRECT HARMONIC BALANCE SOLUTION 
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8.2 BLADE TRANSFER MATRIX RELATIONS 
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BLADE TRANSFER MATRIX RELATIONSHIPS 
“iVshlp^ three important transfer matru; 

the relation of the stat^vMilIbl« an^nfin^m t< l j ai , iable!i in a contmuou! span of the blade (e g 
variables an infiniesimal ‘ ha f tha bl ><“ «P to the s a* 

example). * r,8llt of the outermost hinge, lead-lag in following 

2 ' tt'reult IT T iab ' eS ~ a ouch as a hinge 

with those to the right* theU -^h tge) VariableS t0 ‘ he ° f tha ■-**« “£e 

3 ' ^isplacemem'and'rotat^i st^^rariables at^h^t^ihf tip* * “>* »in g e and the 

formulated eithe^M the sol"^ TIT meth ° d S °'“‘' 0n techni 1"e can be 

direct solution for the harmonic response. ^ ° f ' steady ei l ull ‘ briura Position or as a 


118 



BLADE ELEMENT TRANSFER MATRIX MODEL 


• BLADE STRUCTURAL MODEL IS FULLY COUPLED IN FLAP, 
LEAD-LAG, AND TORSION 


{AZ} t = [Mi(Zi)\{AZ}i+i + {M 2 {Z l , others)} 

where 

{AZ}J = {A M f , A0 f , A F z , AZ, AM t , A0 t , AF y , AY, A M L , A 6 L , A F x , A x } t 

AND {Z t j REPRESENTS THE PREVIOUS ITERATION VALUE 
OF THE STATE VECTOR IN THE STEADY STATE 
SOLUTION AND THE CONVERGED STEADY STATE 
RESPONSE IN THE DYNAMIC PERTURBATION 
SOLUTION 


1 
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BLADE TRANSFER MATRIX RELATIONSHIPS 



HO HI H2 

FLAP TORSION LEAD-LAG 

• TRANSFER RELATIONSHIPS 

{AZ}£ 2 = [<*]{AZ} n+ i + M 
{ &Z} R Hl = [U u }{±Z} L H2 + {« n } 
{ AZ} R H0 = [U 0l \{AZ} L Hl + 


( 1 ) 

( 2 ) 


73 

m 


O 


. discontinuity equations 

{A = [H 2 \{&Z} r H2 - {r fl2 }(A »H2) 
{A Z} L Hi = [ml{AZ}g, - {r"‘}(A9Hl) 


/HERE [HI] AND [H2] ARE KINEMATIC ‘ 

ND A0H2 and A0H1 ARE THE LOCAL ANGLE DISCON 


AT THE HINGES 



r 

* ;, a. 

rn 
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BLADE TRANSFER 


MATRIX RELATIONSHIPS (continued) 



• HINGE EQUILIBRIUM EQUATIONS 

[M 3 ]{ &Z*} n+ i + [M4){A0 d } + [M 5 ] = 0 


(6) 


WHERE 


(AZf, +1 = {A0 f ,AZ,A0 r , AY, Ad L ,AX} N+l 


COMBINING EQUATIONS (1) THROUGH (6) RESULTS ^N THE 
FORMULATION FOR THE ITERATIVE ALGORITHM OF TH 
STEADY STATE SOLUTION OR THE DIRECT SOLUTION OF TH 
HARMONIC PERTURBATION MOTION 
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COUPLED ROTOR/ AIRFRAME FORCED RESPONSE 

(continued) 


• SOLUTION FOR COUPLED ROTOR/ AIRFRAME N-l, N, N+l HARMONIC 
PERTURBATION TIP DISPLACEMENTS AND HINGE DISCONTINUITY 
ANGLES (CONTD) 


where 


• EQUATIONS (7) AND (14) RESULT IN 

(. Eoch\{Vtch } + ^ochKVdch) + {*eoch} = {0} 

[ Eoch ] = 1 B OC h\ - h][/2|(Al[CocH] 

(eocff] = [boon 1 - frPlWleocffl - [CRB\ 

\&eoch} = {S n ocH} - [~i\[R][M{*foch) 

• COMBINE EQUATIONS (12) AND (15) TO YIELD THE SOLUTION 


where 


and 


(15) 

(16) 

(17) 

(18) 


{till] J = -^ a ° ch ^^ aoch} 

(19) 

, . . f [Eoch\ [*och\ 

[Aoch1 - [ [Doom] [docn} \ 


(20) 

1 1 {®eoch} \ 

{oaoch} - { {5doch} J 


(21) 


THE REMAINING BLADE STATE VARIABLES ARE OBTAINED BY 
MULTIPLYING THE TRANSFER MATRICES PROCEEDING FROM TIP TO ROOT 
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9.0 COUPLED ROTOR/AIRFRAME FORCED 
RESPONSE FORMULATION 


127 
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9.1 OVERALL APPROACH 
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OVERALL APPROACH 

of coupled .espouse is the 

matrix method. Airframe impedanee is computed“tag nTsTRAN Thelo^ " si "*‘ he 1 transfer 
the hub are matched to yield the proper counlina r k* AN ‘ The forces and displacements at 

of using impedance matching at the hub is that the 6 a ‘ r ^ rame - The advantage 

one blade. It is assumed that all blades of the r ntr> ° * ! mpedance can be obtained considering only 

position. Since the analysis is perfo^d „ r ‘ ,0ading at the same -tauth 

superposition of the harmonfc ^“es ^ reSp ° nSe is ° btained b Y a 
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OVERALL APPROACH 


• DETERMINED USING THE IMPEDANCE MATCHING 
TECHNIQUE 

• CALCULATION OF ROTOR IMPEDANCE 

• DETERMINATION OF AIRFRAME IMPEDANCE 

• MATCHING OF FORCES AND DISPLACEMENTS AT THE 
HUB 

• ADVANTAGE IN THAT ROTOR IMPEDANCE CAN BE 
OBTAINED BY CONSIDERING A SINGLE BLADE 

• ASSUMPTION INHERENT IN USING IMPEDANCE 
MATCHING: 

• ALL BLADES OF THE ROTOR EXPERIENCE THE SAME 
LOADING AT THE SAME AZIMUTH POSITION 

• RESULTANT RESPONSE IS OBTAINED BY A 
SUPERPOSITION OF HARMONIC PERTURBUATION 
RESPONSES 
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ROTOR IMPEDANCE FORMULATION 
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ROTOR IMPEDANCE FORMULATION 

In order to determine the coupled rotor/airframe response the rotor impedance must be obtained. This 
consists of formulating expressions for the blade root displacements and forces in terms of the 
unknown tip displacements and the angular discontinuities at the hinges. 
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ROTOR IMPEDANCE FORMULATION 


• DETERMINATION OF HARMONIC PERTURBATION BLAD 
ROOT DISPLACEMENTS AND FORCES IN TERMS OF 
PERTURBATION UNKNOWN TIP DISPLACEMENTS AND 
UNKNOWN HINGE ANGULAR DISCONTINUITIES 


{von} = [ b on] {Vtn} + [ 6 o . v ] {»?£*} + {&non} 
{Fqn} = {Con}{vtn} + M{ Vdn) + {*fon} 
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9.2.2 AIRFRAME IMPEDANCE 
CALCULATION 
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AIRFRAME IMPEDANCE CALCULATION 


A NASTRAN finite element model is used to determine the airframe hub mobility. Unit vibratory 
forces and moments are applied to the hub at N/REV in fixed coordinates and the resulting forced hub 
response constitutes the elements of the 6x6 complex response matrix [R] (inverse of the impedance). 
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AIRFRAME IMPEDANCE CALCULATION 


• CALCULATION OF AIRFRAME RESPONSE 

{Z F } = [R]{5 f } 


WHERE THE RESPONSE MATRIX [R] IS DETERMINED FROM 
NASTRAN BY APPLYING UNIT VIBRATION FORCES AND 
MOMENTS TO THE ROTOR HUB AT N/REV IN THE FIXED 
COORDINATE SYSTEM 
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9.2.3 RELATIONSHIP OF ROTATING AND 
FIXED SYSTEM DISPLACEMENTS 

AND FORCES 


143 


relationship op rotating ano™ system disp1acements and 

required in the procraso" d«erroi“nin| coup^ 0 to S // em f diSplaC ' mentS “ d forces ‘he blade root are 
the rotating and fixed system., are related b, cSr gtZo^hT' Displac “ s rotations in 
Forces and moments in the rotating and fixed systems J ? P J a *» b,llty requirements at the blade root, 
individual harmonics in the fixed system. This Lulls j^tiT 6 ^ COnsiderin g equilibrium of the 
harmonics of the rotating system related to N/RFV k the ^ expressions which contain N-l, N, N+l 

forces m the fixed system determine the fuselage respons^) 0 ™ 8 “ **“ fiX6<1 System ’ ( since N/REV 
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RELATIONSHIP OF ROTATING AND FIXED 
SYSTEM DISPLACEMENTS AND FORCES 


• DERIVATION OF COUPLING RELATIONSHIPS FOR 

ROTATING AND NON-ROTATING DISPLACEMENTS AND 
FORCES AT THE BLADE ROOT 

DISPLACEMENTS: {f}^ CH } = \i\{Z F ) + {r/£ p } (10) 

FORCES: {S F } = [A }{Fg CH } ( n ) 
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.4 INBOARD HINGE MOMENT 
EQUILIBRIUM EQUATIONS 



INBOARD HINGE MOMENT EQUILIBRIUM EQUATIONS 


M°ment equilibrium equations at the hinges are required for the N-l N and 
equilibrium equations are expressed in terms of the harmonic tip state vector 
rotations and the harmonic discontinuities in the angles across the hinges. 


N+l harmonics. These 
displacements and 
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INBOARD HINGE MOMENT EQUILIBRIUM EQUATIONS 


• MOMENT EQUILIBRIUM EQUATIONS AT HINGES IN 

TERMS OF TIP DISPLACEMENTS AND DISCONTINUITIES 
ACROSS THE HINGES 


[ DochKvtch } + l d 0CH\{r)DCH} + {&doch} = { 0 } ( 12 ) 
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9.2.5 BOUNDARY CONDITION 
FORMULATION USING HARMONIC 
ANALYSIS APPROACH 
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BOUNDARY CONDITION FORMULATION USING HARMONIC ANALYSIS 

APPROACH 


Using the frequency domain-harmonic analysis approach enables a clear distinction and 
straightforward formulation of the boundary conditions required for the solution of the coupled 
rotor/airframe response problem. For the steady equilibrium position solution, the rotor hub is 
considered to be rigid. For an N-bladed rotor the N-l, N, and N+l harmonics are coupled with the 
airframe using impedance matching. The remaining harmonics are assumed to be uncoupled from the 
airframe motion and are solved for separately. 
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BOUNDARY CONDITION FORMULATION USING 
HARMONIC ANALYSIS APPROACH 


• STEADY STATE SOLUTION IS DETERMINED USING AN 
ITERATIVE APPROACH ASSUMING A RIGID ROTOR HUB 


• ( N_1 )» N, (N+l) HARMONIC (N= NUMBER OF ROTOR 
BLADES) SOLUTION COUPLED WITH AIRFRAME 

FROM 

• OTHER HARMONICS UNCOUPLED^AIRFR AME MOTION 
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9.3 COUPLED ROTOR/AIRFRAME FORCED 

RESPONSE 
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COUPLED ROTOR/ AIRFRAME FORCED RESPONSE 

The solution for the coupled rotor/airframe forced response is obtained in a four-step procedure: 

1 J he N_1 ’ N ’ and N+1 harmonic rotating root displacements are expressed in terns of the 

K T nd f ° rCe rotatin 8 _to ' nonr otating transformation matrices, [7J and [A], respectively, 
he NASTRAN response matrix [Rj, and the N-l, N, and N+1 harmonic rotating blade root forces 
(Equation ( 13 )). 

2. The expression for the N-l, N and N+1 harmonic rotating blade root forces (Equation (8)) which 
is in terms of the harmonic blade tip displacements and rotations and hinge angular 
discontinuities, is substituted into the expression for the N-l, N, and N+1 harmonic rotating root 
displacements m 1) above. This results in the N-l, N, and N+1 rotating root displacements 
expressed in terms of the N-l, N, and N+1 harmonic tip displacements and hinge angular 
discontinuities (Equations ( 14 )). 
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COUPLED ROTOR/ AIRFRAME FORCED RESPONSE 


• SOLUTION FOR COUPLED ROTOR/ AIRFRAME N-l, N, N+i 
HARMONIC PERTURBATION TIP DISPLACEMENTS AND 
HINGE DISCONTINUITY ANGLES 

• COMBINING EQUATIONS (9), (10), and (11) 


inScH) = iimwvScH} + iti p } (i3) 


• SUBSTITUTE FOR THE ROTATING VECTOR {Fq CH } FROM 
EQUATION (8) 


{Voch} - [l\[ R \m[CocH]{VTCH} + bPIMI Coch]{t1tch } (14) 

+ fr][tf]['M{ < *FOC7if} + {Vd P } 
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COOPUD KOTOK/AlRPItAMB porcbo pesp„„ c 

are equated. + l>armonic tlp displacements and hinge discontinuities 

equations a ^ 7* N+1 mom ™‘ ^librium 

the N - *- d "« 

In order to determine the remaining N-l N 

transfer matrices and discontinue relation' ace — * 
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COUPLED ROTOR/ AIRFRAME FORCED RESPONSE 

(continued) 


• SOLUTION FOR COUPLED ROTOR/ AIRFRAME N-l, N, N+l HARMONIC 
PERTURBATION TIP DISPLACEMENTS AND HINGE DISCONTINUITY 
ANGLES (CONTD) 

• EQUATIONS (7) AND (14) RESULT IN 

[Eoch\{^tch) + l e oc«l{»»£cjj} + {&eoch} = {0} 

w h ere 

[ E 0 ch\ = ( Boch ] ~ ] 

[coch] = [ 1>och \ - - [CRB] 

[ oleoch | = {«r jOCh) - (7|[-R][^]{“foch} 


• COMBINE EQUATIONS (12) AND (15) TO YIELD THE SOLUTION 



{ it;) } - 

(19) 

where 

r a I _ [ l-^OCffl [coch\ 

[A ° ch1 - [ [Doch\ [doca ]J 


(20) 

and 

- { tel } 


(21) 

THE REMAINING BLADE 

STATE VARIABLES ARE OBTAINED BY 



MULTIPLYING THE TRANSFER MATRICES PROCEEDING FROM TIP TO ROOT 
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ANALYSIS AND CORRELATION PLAN 


The analysis and correlation plan is designed to follow f„ u , steps 

(i) ~ - - — —ary 

(li) m ° del * ‘—ting a forced 

<iH) Bn “ e eleme "‘ SlrUCtUral m ° dcl b >- “" d -ting several checks on the fuselage 

' mod “"Ut SXic^lt^'sitgle mothiufi program^ 

modeling P ’.rrora coiljatton rtS** *\ a Sp * ciflc fli * ht co " dili °" (1M knots) to identify 

Satisfactory correlation will be obtained at IK knoTs p^rVo “ alylu'a t “oter B^rcond” 015 ’' 
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ANALYSIS AND CORRELATION PLAN 


• DERIVE BLADE STRUCTURAL MODEL 

• VERIFY BLADE STRUCTURAL MODEL THROUGH MODAL 
SURVEY 

• VERIFY FUSELAGE FINITE ELEMENT MODEL 

• DEFINE AND INTEGRATE AERODYNAMIC MODEL IN 
RACAP 

• CORRELATE RACAP AND TEST DATA FOR ONE FLIGHT 
CONDITION (114 KNOTS) 

• PERFORMANCE ANALYSIS AT OTHER FLIGHT 
CONDITIONS 
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10.1 AH-1G BLADE STRUCTURAL 

MODELING 
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AH-1G BLADE STRUCTURAL MODELING 


In the frequency domain, where a single blade is modeled to simulate an entire rotors’ response, the 
teetering rotor is modeled using different boundary conditions for even and odd harmonics. These are, 
cantilevered for even, and pinned for odd. The model includes underslinging and precone effects. The 
lag hinge is absent, while the control system is modeled by a single spring. The model is limited by the 
assumptions listed below. 

• Fuselage vibrations are uninfluenced by control system/aerodynamic interaction. 

• The fuselage FEM is linear. 

• The elastomeric “soft” transmission support is independent of frequency. 

• The higher frequency vibrations are not calculated (> 2 rev). 
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AH-1G BLADE STRUCTURAL MODELING 


• MODEL EVEN AND ODD HARMONIC BOUNDARY 
CONDITIONS INDIVIDUALLY 

• INCLUDE UNDERSLINGING AND PRECONE EFFECTS 

• LAG HINGE - ABSENT 

• CONTROL SYSTEM MODELED BY SINGLE SPRING 
ASSUMPTIONS 

• FUSELAGE VIBRATION THROUGH CONTROL SYSTEM 
AND AERODYNAMIC INTERACTION NEGLECTED 

• NON LINEARITIES OF FUSELAGE MODEL NEGLECTED 

• LIMITED ACCESS TO PROPERTIES OF “SOFT” 
TRANSMISSION SUPPORT 

• FUSELAGE VIBRATION PREDICTED AT 2 /REV (4/REV, 
6/REV NEGLECTED) 
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AH-1G AERODYNAMICS MODEL 



AH-1G AERODYNAMICS MODEL 


An isolated rotor trim program was used to define the rigid blade airloads using (a) Glauert inflow 
and (b) a free wake inflow model. Both analyses used a blade element analysis (lifting line model) 
and experimentally measured 2-D airfoil data. Airloads correlations with flight test data were 
done with both inflow models. The free wake model was adopted for the remaining analysis 
because it yielded better correlation with measured airloads. 
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AH-1G AERODYNAMICS MODEL 


• ISOLATED ROTOR TRIM 

• GLAUERT/FREE WAKE INFLOW MODEL 

• RIGID BLADE ANALYSIS 

• BLADE ELEMENT ANALYSIS 

• EXPERIMENTALLY MEASURED STEADY 2D AIRFOIL 
DATA 
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10.3 


AH-1G AIRFRAME MODEL 
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AH-1G AIRFRAME MODEL 


Several checks were conducted on the fuselage finite element model to ensure its validity. These 
were, (i) the multi-level strain energy check, (ii) kinetic energy check, (iii) connectivity check. In 
the multi-level strain energy check, the strain energy of the free-free model are calculated and 
used to identify errors in modelling in the “G”, “N”, “F” set levels. In the kinetic energy check, 
global and local structural modes can be identified with computational ease and speed. The 
connectivity check identifies areas of singularity in the model due to misaligned or incomplete 
connectivities in the model. 

The hub impedance matrix (at 2/Rev) is used in arriving at coupled rotor-fuselage hub loads. This 
matrix was calculated using the airframe NASTRAN model, employing a modal response solution 
sequence. Transfer matrices (at 2/Rev) relating unit hub loads to vibrations at locations of 
interest on the fuselage are generated using the same airframe model and modal response solution. 
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AH-1G AIRFRAME MODEL 


• CHECKS CONDUCTED ON AH-1G FEM: 

• MULTI-LEVEL STRAIN ENERGY 

• KINETIC ENERGY 

• CONNECTIVITY 

• HUB IMPEDANCE MATRIX GENERATED AT 2 REV 
(MODAL RESPONSE SOLUTION) 

• FUSELAGE VIBRATION RESPONSE MATRICES 
GENERATED AT 2P (MODAL RESPONSE SOLUTION) 
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10.4 COUPLED ROTOR-FUSELAGE 

ANALYSIS 
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COUPLED ROTOR-FUSELAGE ANALYSIS 


The coupled rotor-fuselage procedure was described in detail earlier for a generic helicopter. For 
the AH-1G, the matched impedance method is applied to the 1, 2, and 3 harmonics of rotor RPM. 
All other harmonics up to 9 are modelled with a rigid hub (i.e. uncoupled from the fuselage 
motions). This assumption makes use of the fact that the rotor acts as a frequency filter and 
transmits only loads at integral multiples of the blade passage frequency to the fuselage. The hub 
loads obtained in the rotating system are transformed to the fixed system at blade passage 
frequency and are then used with the fuselage transfer matrices to compute the fuselage vibration 
response at locations of interest. 
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COUPLED ROTOR-FUSELAGE ANALYSIS 


• AIRFRAME HUB IMPEDANCE FROM FEM MATCHED 
WITH ROTOR HUB IMPEDANCE 

• 1, 2, 3, HARMONICS ANALYSIS PERFORMED WITH 
FUSELAGE COUPLING 

• ALL OTHER HARMONICS (UP TO 9) ARE UNCOUPLED 
FROM FUSELAGE 

• 2P HUB LOADS (IN FIXED SYSTEM) APPLIED TO 
FUSELAGE RESPONSE MATRICES FOR VIBRATION 
LEVEL CORRELATION 


179 


180 



DISCUSSION OF RESULTS 



DISCUSSION OF RESULTS 


Fuselage acceleration levels (lateral and vertical) for all flight conditions (steady forward flight, 67 
to 142 knots) compare favorably. Some general trends are identified below. 

(a) Vertical vibration prediction is generally better correlated than lateral vibration predictions. 
This trend is consistent with the results of the study of the finite element model and the test data 
(Appendix 1). Those results indicate a factor of six discrepancy between calculated (derived from 
measured mast top accelerations at 85 knots) and measured lateral acceleration values at the nose. 
In that study, the vertical vibration levels are much more consistent than the lateral. RACAP 
predictions and measured vibration levels show a discrepancy of about a factor of 6.5 in the lateral 
hub vibration prediction at 85 knots. This discrepancy feeds into the vibration prediction at any 
given location on the ship in the lateral direction. However, the problem appears to lie in the 
finite-element formulation rather than the load prediction methodology. In Appendix 1, the 
fuselage finite element model effectively has a higher impe ance at 2P than is realistic. This 
manifests itself in two ways: (i) the root boundary conditions are overly stiff and therefore the hub 
loads are magnified; and (ii) for a given hub load the fuselage accelerations are underestimated. 

Of these, the magnification of coupled rotor-fuselage hub loads is probably more pronounced due 
to the strong dependence of the response on the fuselage impedance. Insofar as this behavior is 
more pronounced in the lateral direction, correlation of fuselage acceleration in that direction is 
expected to be poorer than the vertical, although the latter is also affected. 

(b) Lower speed flight conditions generally show better correlation than higher speed conditions at 
the same location. 

The nonlinearity of the elastomeric mounts is not adequately characterized in the 
longitudinal/lateral directions. The nonlinear stiffness and damping of the mounts varies with 
'oad magnitude; this effect is not ncluded in the transfer matrix formulation. At higher flight 
speeds, the hub impedance r "' -t: defining the root boundary conditions would also be changed 

by f he inclusion of this efTcr’’"' Hence the inconsistency. 

(c) I irward locations on ip indicate better correlation than aft locations. 

A major inflm nee on aft k ution vibrations is the 2P main rotor wake impinging on the aft 
fuselage. Thes, k-ds are not a part of RACAP formulation. Since these loads are not calculated 
in the present k rmulation, the influence of this effect on the vibratory response is not accounted 
for. In addition, 2P control system link loads have an influence on the fuselage vibrations that are 
not part o, this analysis. These vibration predictions may have been improved by the inclusion of 
these effects. 
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DISCUSSION OF RESULTS 


• VERTICAL FUSELAGE VIBRATION LEVELS AFFORD 
GOOD CORRELATION WITH MEASURED DATA IN 
FORWARD PART OF SHIP 

• AFT VERTICAL VIBRATION CORRELATION NOT AS 
GOOD DUE TO MISSING EMPENNAGE ROTOR WAKE 
EXCITATION 

• LATERAL VIBRATION RESPONSE CORRELATION NOT AS 
GOOD AS VERTICAL VIBRATION LEVEL CORRELATION 

• FUSELAGE FEM NEEDS REFINEMENT - CHECK NOT 
SATISFACTORILY CONSISTENT (REF. APPENDIX 1) 
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RESULTS 


Nose 2P (Lateral and Vertical) Acceleration level against flight speed 
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RESULTS 


Pilot Seat 


2P (Lateral and Vertical) Acceleration level against 
flight speed 
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RESULTS 


Pilot Seat 2P (Lateral and Vertical) Acceleration level against 
flight speed 






RESULTS 


Hub 
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RESULTS 


Hub 2P (Lateral and Vertical) Acceleration level against 
flight speed 
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RESULTS 


CG 2P (Lateral and Vertical) Acceleration level against 
flight speed 
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RESULTS 


CG 2P (Lateral and Vertical) Acceleration level against 
flight speed 



RESULTS 


Gear box 


2P (Lateral and Vertical) Acceleration level against 
flight speed 


no lateral test data available 
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RESULTS 

Gear box 2P (Lateral and Vertical) Acceleration level against 
flight speed 

no lateral test data available 
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Engine 
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flight speed 
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RESULTS 


Engine 2P (Lateral and Vertical) Acceleration level against 
flight speed 




RESULTS 


R-Wing 2P (Lateral and Vertical) Acceleration level against 
flight speed 

no test data available 
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RESULTS 


L-Wing 


2P (Lateral and Vertical) Acceleration level against 
flight speed 


no test data available 
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RESULTS 


T-B Junction 
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RESULTS 


Tail 2P (Lateral and Vertical) Acceleration level against 
night speed 

no lateral test data available 
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RESULTS 


Tail 2P (Lateral and Vertical) Acceleration level against 
flight speed 


no lateral test data available 
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RESULTS 


Elevator 


2P (Lateral and Vertical) Acceleration level against 
flight speed 
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RESULTS 


Elevator 2P (Lateral and Vertical) Acceleration level against 
flight speed 



RESULTS 


Fin 2P (Lateral and Vertical) Acceleration level against 
flight speed 
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CONCLUDING REMARKS 


A description and summary of MDHC’s method for coupled rotor/airframe analysis, RACAP has 
been presented m the above discussion. RACAP is capable of accurately modeling all rotor 

for a ^v^T hr U ^ tee ^' ng r ° t0rS> t0 determine rotor and Pelage forced response loads 
mod,, I f COnd,t,on \ The ana| y sis has formulated to facilitate the use of different 

modules for inflow, rotor trim, unsteady aerodynamic airload calculations, and structural model 
root end boundary conditions. The modular approach allows the analyst to easily update the 

LT, am t investigate alternate analytical approaches to a particular part of the overall 

aeroelastic response problem. 

A description of the results obtained using RACAP to model the AH-1G has also been presented. 

wake infl ^ ^ ^ trends ’ there are some ^ uestions regarding the input data. Free 

ake inflow generally indicated better correlation. Inclusion of wake-fuselage interaction in he aft 
regions may improve correlation. ne ait 


In general, fairly good correlation is obtained for all flight conditions. 
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CONCLUDING REMARKS 


• MDHC’S COUPLED ROTOR/ AIRFRAME ANALYSIS, 

RACAP IS A COMPREHENSIVE AEROELASTIC PROGRAM 
THAT CAN ADEQUATELY MODEL ALL ROTOR 

C ONFIGUR ATION S 

• RACAP HAS BEEN DEVELOPED WITH THE FLEXIBILITY 
TO INCORPORATE DIFFERENT STRUCTURAL AND 
AERODYNAMIC MODULES 

• AH-1G MODELING SHOWS CONSISTENT AND FAIRLY 
ACCURATE FUSELAGE VIBRATION LEVEL 
CORRELATION 
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In order to evaluate the validity of the AH-1G Finite Element Model from a coupled rotor-fuselage 
vibration prediction standpoint, an analytical test was designed and performed. The test 
consisted of using the measured hub vibrations, and computing from these an estimate of the hub 
loads required to produce such vibrations. This relationship is 


Fx 

F y 

Fz 



(i) 


In equation 1, the left hand side is an estimate of the hub loads required to produce the vibration 
level measured on the ship, the matrix R„ H is a matrix calculated from the fuselage finite element 
model, relating unit hub forces to hub vibrations, and the vector on the right hand side is the 
measured hub accelerations. By applying these hub forces to a transfer matrix relating hub forces 
to fuselage vibrations at any location of interest, L, one can compute the vibratory response at L. 
If the fuselage finite element model were exact, this vibration level at L would be equal to the 
measured quantity. This test was performed at all the locations at which test data was available, 
and did not provide good correlation at any of these locations. In particular, the lateral vibration 
levels are most in error, differing by up to 80 percent in magnitude. The vertical vibration 
magnitudes are significantly better, but still not adequate. These results lead to the conclusion 
that the fidelity of the NASTRAN model is less than adequate, or that there are additional 
factors affecting vibrations, not included in the model. The RACAP predicted hub loads and 
vibrations depend on the hub impedance as well as the fuselage transfer matrices, both of which 
are a function of the finite element model. 

In conclusion, there are reasons to suspect the validity of the NASTRAN model, and these 
provide some of the reasons for errors in the predicted fuselage vibratory response. 
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